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Abstract 

The interaction of a small vortex ring with the free surface of a perfect fluid 
is considered. In the frame of the point ring approximation the asymptotic 
expression for the Fourier-components of radiated surface waves is obtained in 
the case when the vortex ring comes from infinity and has both horizontal and 
vertical components of the velocity. The non-conservative corrections to the 
equations of motion of the ring, due to Cherenkov radiation, are derived. 

1 Introduction 

The study of interaction between vortex structures in a fluid and the free surface is 
important both from practical and theoretical points of view. In general, a detailed 
investigation of this problem is very hard. Even the theories of potential surface 
waves and the dynamics of vortices in an infinite space taken separately still have 
a lot of unsolved fundamental problems on their own. Only the consideration of 
significantly simplified models can help us to understand the processes which take 
place in the combined system. 

In many cases it is possible to neglect the compressibility of the fluid as well as 
the energy dissipation. Therefore the model of ideal homogeneous incompressible 
fluid is very useful for hydrodynamics. Because of the conservative nature of this 
model the application of the well developed apparatus of Hamiltonian dynamics be- 
comes possible |T[ [[|. An example of effective use of the Hamiltonian formalism in 
hydrodynamics is the introduction of canonical variables for investigations of poten- 
tial flows of perfect fluids with a free boundary. V.E.Zakharov showed at the end 
of the sixties that the surface shape z = T](x,y,t) and the value of the velocity 
potential ip(x,y,t) on the surface can be considered as generalized coordinate and 
momentum, respectively. 
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It is important to note that a variational formulation of Hamiltonian dynamics in 
many cases allows to obtain good finite-dimensional approximations which reflect the 
main features of the behavior of the original system. There are several possibilities 
for a parameterization of non-potential flows of perfect fluid by some variables with 
dynamics determined by a variational principle. All of them are based on the con- 
servation of the topological characteristics of vortex lines in ideal fluid flows which 
follows from the freezing-in of the vorticity field ft(r, t) = curlv(r, t). In particular, 
this is the representation of the vorticity by Clebsch canonical variables A and fi || 

B 

n(r,*) = [VA x V/x] 

However, the Clebsch representation can only describe flows with a trivial topology 
(see, e.g., ||). It cannot describe flows with linked vortex lines. Besides, the variables 
A and \x are not suitable for the study of localized vortex structures like vortex 
filaments. In such cases it is more convenient to use the parameterization of vorticity 
in terms of vortex lines and consider the motion of these lines , [0] , even if the global 
definition of canonically conjugated variables is impossible due to topological reasons. 

This approach is used in the present article to describe the interaction of deep (or 
small) vortex rings of almost ideal shape in the perfect fluid with the free surface. In 
the case under consideration the main interaction of the vortex rings with the surface 
can be described as the dipole-dipole interaction between "point" vortex rings and 
their "images". Moving rings interact with the surface waves, leading to radiation 
due to the Cherenkov effect. Deep rings disturb the surface weakly, so the influence 
of the surface can be taken into account as some small corrections in the equations 
of motion for the parameters of the rings. 

In Sec. 2 we discuss briefly general properties of vortex line dynamics, which fol- 
low from the freezing-in of the vorticity field. In Sec. 3 possible simplifications of 
the model are made and the point ring approximation is introduced. In Sec. 4 the 
interaction of the ring with its image is considered. In Sec. 5 we calculate the Fourier- 
components of Cherenkov surface waves radiated by a moving vortex ring and deter- 
mine the non-conservative corrections caused by the interaction with the surface for 
the vortex ring equations of motion. 

2 Vortex lines motion in perfect fluid 

It is a well known fact that the freezing-in of the vorticity lines follows from the Euler 
equation for ideal fluid motion 

ft t = curl [v x ft], v = curl -1 ft 

Vortex lines are transported by the flow |L| , |4[] , |8|] . They do not appear or disappear, 
neither they intersect one another in the process of motion. This property of perfect 
fluid flows is general for all Hamiltonian systems of the hydrodynamic type. For 
simplicity, let us consider temporally the incompressible fluid without free surface in 
infinite space. The dynamics of the system is specified by a basic Lagrangian L[v], 
which is a functional of the solenoidal velocity field. The relations between the ve- 
locity v, the generalized vorticity ft, the basic Lagrangian L[v] and the Hamiltonian 
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H[Q] are the following [|Q 



n = curl f ^ ] v = v [n] 



5v / 



(2) 

v=v[n] v ; 



v = curl (3) 

and the equation of motion for the generalized vorticity is 

n t = curl [curl (6H/6Q) x fi] (4) 

This equation corresponds to the transport of frozen-in vortex lines by the velocity 
field. In this process all topological invariants |T0 of the vorticity field are conserved. 
The conservation of the topology can be expressed by the following relation 

n ( r,t) = fnr- tt(M))W.)v.WM)fc^™| Mi)| (5) 

where the mapping R(a, t) describes the deformation of lines of some initial solenoidal 
field fio(r). Here a(r, t) is the inverse mapping with respect to R(a, t). The direction 
of the vector b 

b(a,t) = (n (a)V a )R(a,t) (6) 

coincides with the direction of the vorticity field at the point R(a, t). The equation 
of motion for the mapping R(a, t) can be obtained with the help of the relation 

fi t (r, t) = curl r J 5(r - R(a, t)) [R t (a, t) x b(a, t)]da, (7) 
1 For the ordinary ideal hydrodynamics in infinite space the basic Lagrangian is 

f v 2 

LsuierM = / — dr fi = curlv 

The Hamiltonian in this case coincides with the kinetic energy of the fluid and in terms of the 
vorticity field it reads 

HEuierM = -1/2 / nA-'ndr = A- // ^MlilM dndra 
J 8vr J J |ri-r 2 | 

where A -1 is the inverse Laplace operator. 

Another example is the basic Lagrangian of Electron Magneto-hydrodynamics which takes into 
account the magnetic field created by the current of electron fluid through the motionless ion fluid. 

LemhdM = - /v(l - A~Vdr => fTl = curl(l - A _1 )v 

n E MH D [n} = iy n(i - A)- 1 ^ = ^JJ ^~^r n ( ri ) ' n ( r 2)*idr 2 

The second example shows that the relation between the velocity and the vorticity can be more 
complex than in usual hydrodynamics. 
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which immediately follows from Eq.(|5]). The substitution of Eq.([7|) into the equation 
of motion (||) gives [|Tl"11 



curl. 



'b(a,t) x [Rt(a,t) - v(R,t)] > 
det||«9R/<9a|| 







One can solve this equation by eliminating the curl r operator. Using the general 
relationship between variational derivatives of some functional F[Q] 



b x curl 



5F 



sn(R) 



5F 



5R(a) 



it is possible to represent the equation of motion for R(a, t) as follows 



[(n (a)V a )R(a) x R t (a)] 



6H\n\R] 



5R(a) 



n 



(9) 



It is not difficult to check now that the dynamics of the vorticity field with topological 
properties defined by Oo in the infinite space is equivalent to the requirement of an 
extremum of the action (5S = 5 J Cn dt = 0) where the Lagrangian is [[5J 



Cn = \j ([R t (a) x R(a)] ■ (fi (a)V a )R(a))da - H[Q[R]]. 



(10) 



In the simplest case, when all vortex lines are closed it is possible to choose new 
curvilinear coordinates u±, u 2 , £ in a-space such that Eq.(^) can be written in a simple 
form 

n(r, t) = J d 2 u f 5(r - R(u, £, t))R^. (11) 

Here v is the label of a line lying on a fixed two-dimensional manifold A/", and £ 
is some parameter along the line. It is clear that there is a gauge freedom in the 
definition of v and £. This freedom is connected with the possibility of changing the 
longitudinal parameter £ = v, t) and also with the relabeling of v 



v 



:i2) 



Now we again consider the ordinary perfect fluid with a free surface. To describe 
the flow entirely it is sufficient to specify the vorticity field fi(r, t) and the motion of 
the free surface. Thus, we can use the shape R(^, £, t) of the vortex lines as a new 
dynamic object instead of f2(r, t). It is important to note that in the presence of the 
free surface the equations of motion for R(u, £, t) follow from a variational principle 



as in the case of infinite space. It has been shown |L2] that the Lagrangian for a 
perfect fluid, with vortices in its bulk and with a free surface, can be written in the 
form 

£ = ||d 2 ^([R t xR].R^ + J*ntdr ± -H\R,*,7j\. (13) 

The functions ^(rj_, t) and r)(r±, t) are the surface degrees of freedom for the system. 
\& is the boundary value of total velocity potential, which includes the part from vor- 
tices inside the fluid, and 77 is the deviation of the surface from the horizontal plane. 
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This formulation supposes that vortex lines do not intersect the surface anywhere. 
In the present paper only this case is considered. 

The Hamiltonian TC in Eq. ([13]) is nothing else than the total energy of the system 
expressed in terms of [R, rf\. 

Variation with respect to R(^, £, t) of the action defined by the Lagrangian ( |13|) 
gives the equation of motion for vortex lines in the form 

[R{ x R(] = ffl[npi,«,i] . (14) 

This equation determines only the transversal component of R 4 which coincides with 
the transversal component of the actual solenoidal velocity field. The possibility of 
solving Eq. flli]) with respect to the time derivative Rt is closely connected with the 
special gauge invariant nature of the TC[R] dependence which results in 

m R -n 

The tangential component of Rj with respect to vorticity direction can be taken 
arbitrary. This property is in accordance with the longitudinal gauge freedom. The 
vorticity dynamics does not depend on the choice of the tangential component. 

Generally speaking, only the local introduction of canonical variables for curve 
dynamics is possible. For instance, a piece of the curve can be parameterized by one 
of the three of Cartesian coordinates 

R={X{z,t),Y{z,t),z) 

In this case the functions X(z,t) and Y(z,t) are canonically conjugated variables. 
Another example is the parameterization in cylindrical coordinates, where variables 
Z(9,t) and (1/2)R 2 (9, t) are canonically conjugated. 

Curves with complicated topological properties need a general gauge free descrip- 
tion by means of a parameter £. 

It should be mentioned for clarity that the conservation of all vortex tube volumes, 
reflecting the incompressibility of the fluid, is not the constraint in this formalism. It 
is a consequence of the symmetry of the Lagrangian (|13"D with respect to the relabel- 
ing (|12D v — > v 0. Volume conservation follows from that symmetry in accordance 
with Noether's theorem. To prove this statement, we should consider such subset of 
relabelings which forms a one-parameter group of transformations of the dynamical 
variables. For small values of the group parameter, r, the transformations are de- 
termined by a function of two variables T(u\, v 2 ) (with zero value on the boundary 
dM) so that 

/ ffr fix \ 

RKf2,^^i,^^Rk-^+0(T 2 ), u 2 + r— + 0(r 2 ) 1 e) (15) 



Due to Noether's theorem, the following quantity is an integral of motion [|1S 
8C <9R^ 

T=0 
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After simple integrations in parts the last expression takes the form 

I T = f A/r(i/i,i/ 2 )([R 1 xR 2 ].R { K= / T{v l ,v 2 )V{v l ,v 2 ,t)d 2 v (16) 

where V(z/ 1; u 2 , t)d 2 v is the volume of an infinitely thin vortex tube with cross-section 
d 2 v. It is obvious that actually the function V doesn't depend on time t because the 
function T(z/ 1; u 2 ) is arbitrary Q, 

3 Point ring approximation 

In general case an analysis of the dynamics defined by the Lagrangian ([ID is too 
much complicated. We do not even have the exact expression for the Hamiltonian 
7i\R, ty,r)] because it needs the explicit knowledge of the solution of the Laplace 
equation with a boundary value assigned on a non-flat surface. Another reason is 
the very high nonlinearity of the problem. 

In this paper we consider some limits where it is possible to simplify the system 
significantly. Namely, we will suppose that the vorticity is concentrated in several 
very thin vortex rings of almost ideal shape. For a solitary ring the perfect shape is 
stable for a wide range of vorticity distributions through the cross-section. This shape 
provides an extremum of the energy for given values of the volumes of vortex tubes 
and for a fixed momentum of the ring. As already mentioned, volume conservation 
follows from Noether's theorem. Therefore some of these quantities (those of which 
are produced by the subset of commuting transformations) can be considered as 
canonical momenta. Corresponding cyclical coordinates describe the relabeling ([12]) 
of the line markers, which doesn't change the vorticity field. Actually these degrees 
of freedom take into account a rotation around the central line of the tube. This 
line represents the mean shape of the ring and we are interested in how it behaves in 
time. For our analysis we don't need the explicit values of cyclical coordinates, but 
only the conserved volumes as parameters in the Lagrangian. 

A possible situation is when a typical time of the interaction with the surface and 
with other rings is much larger then the largest period of oscillations corresponding 
to deviations of the ring shape from perfect one. Under this condition, excitations 
of all (non-cyclical) internal degrees of freedom are small during all the time, and a 
variational anzats completely disregarding them reflects the behavior of the system 
adequately. The circulations 

r„ = / d 2 u 

of the velocity for each ring don't depend on time. A perfect ring is described by the 
coordinate R n of the center and by the vector P n = r n S n , where S n is an oriented 
area of the ring. We use in this work the Cartesian system of coordinates (x, y, z), 
so that the vertical coordinate is z, and the unperturbed surface is at z = 0. The 

2 If vortex lines are not closed but form a family of enclosed tori then the relabeling freedom is 
less rich. In that case one can obtain by the similar way the conservation laws for volumes inside 
closed vortex surfaces. Noether's theorem gives integrals of motion which depend on an arbitrary 
function of one variable S(Q, where £ is the label of the tori. 
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corresponding components of the vectors R n and P n are 

R71 Y ny Znj : P n {^Pxni Pyni Pzn) 

It is easy to verify that the vectors P n are canonically conjugated momenta for the 
coordinates R n . To verify that we can parameterize the shape of each vortex line in 
the following manner 

M 

R(f,t)= £ r m (t)e im «, r_ m = f m (17) 

m=—M 

Here r m (t) are complex vectors. Substituting this into the first term of the La- 
grangian flT3|) gives 

- <f(\Rt x R]R c )d£ = 2 7 rzf ([r_ 1 x n] + 2[r_ 2 x r 2 ] + . . .) + 

+ + 27rir_ 1 [r_i x r 2 ] - 27rzri[ri x r_ 2 ] + • • • (18) 

If we neglect the internal degrees of freedom which describe deviations of the ring 
from the ideal shape 

(r_x) 2 = (n) 2 = 0, r 2 = r_ 2 = 0, 

then the previous statement about canonically conjugated variables becomes obvious: 

R n = r 0n , P„ = 2nT n ■ z[r_ ln x r ln ] (19) 

Such an approximation is valid only in the limit when sizes of rings are small in 
comparison with the distances to the surface and the distances between different 
rings 

4/^<|Z n |,|R n -R i |, l^n. (20) 

V 1 n 

These conditions are necessary for ensuring that the excitations of all internal degrees 
of freedom are small. Obviously, this is not true when a ring approaches the surface. 
In that case one should take into account also the internal degrees of freedom for the 
vortex lines. 

The inequalities ( p0|) also imply that vortex rings in the limit under consideration 
are similar to point magnet dipoles. This analogy is useful for calculation of the 
Hamiltonian for interacting rings. In the main approximation we may restrict the 
analysis by taking into account the dipole-dipole interaction only. 

It should be mentioned that in some papers (see e.g. |14| and references in that 
book) the discrete variables identical to R n and P n are derived in a different way 
and referred as the vortex magnetization variables. 

In the expression for the Hamiltonian, several simplifications can be made. Let 
us recall that for each moment of time it is possible to decompose the velocity field 
into two components 

v = V o + V0. (21) 
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Here the field V satisfies the following conditions 



(V-V o ) = 0, curlV = n, (n-V )| 2= „ = 0. 

The boundary value of the surface wave potential 0(r) is ^(rj_). In accordance with 
these conditions the kinetic energy is decomposed into two parts and the Hamiltonian 
of the fluid takes the form 

n = \ J V 2 dh + ^( V0 -dS) + 9 -J V 2 dr ± (22) 

The last term in this expression is the potential energy of the fluid in the gravitational 
field. If all vortex rings are far away from the surface then its deviation from the 
horizontal plane is small 

|Vt7|«1, \v\^\Z n \ (23) 

Therefore in the main approximation the energy of dipoles interaction with the sur- 
face can be described with the help of so called "images". The images are vortex 
rings with parameters 

R n (^nj ^ti) ^n)i P n (Pxm Pyni Pzn) (^4) 

The kinetic energy for the system of point rings and their images is the sum of the 
self-energies of rings and the dipole-dipole interaction between them. The expres- 
sion for the kinetic energy of small amplitude surface waves employs the operator 
k which multiplies Fourier-components of a function by the absolute value A; of a 
two-dimensional wave vector k. So the real Hamiltonian H is approximately equal 
to the simplified Hamiltonian 7i 

n~H = l [ (#V> + 9V 2 )dr± + Y. £ n(Pn) + 

J n 

1 ^ 3((R n - R|) • P n )((R n - R*) • P*) - |R n - Rf | 2 (P n • Pf) 

+ ^WF (25) 

With the logarithmic accuracy the self-energy of a thin vortex ring is given by 
the expression 

^"T^A^r) (26) 

where the small constant A n is proportional to the conserved volume of the vortex 
tube forming the ring. This expression can easily be derived if we take into account 
that the main contribution to the energy is from the vicinity of the tube where 
the velocity field is approximately the same as near a straight vortex tube. The 
logarithmic integral should then be taken between the limits from the thickness of 
the tube to the radius of the ring. 



3((R n — Ri) • Pn)((R« 


- Ri) P 


*)-! 


Rn — R« 


2 (Pn • P,) , 




Rn — R« 


5 1 
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In the relation ^ = $o+V ; the potential $o is approximately equal to the potential 
created on the flat surface by the dipoles and their images 

%(r x ) « $(r ± ) = -i- £ (P ; r (r _ ± R ^ (27) 

i/l n |IJ_ -tv-n| 

In this way we arrive at the following simplified system describing the interaction 
of point vortex rings with the free surface 

C = £ R n P n + [r](iP + $)d 2 r ± - W[{Rn, P„}, V, V>] (28) 

It should be noted that due to the condition (|20|) the maximum value of the 
velocity Vo on the surface is much less then the typical velocities of the vortex rings 

p r 3 / 2 
" « 



Z 3 pV2 

Therefore the term V 2 /2 in the Bernoulli equation 

+ Vo + 9V + small corrections = 

is small in comparison with the term ty t - The Lagrangian (|28|) is in accordance with 
this fact because it does not take into account terms like (1/2) / V 2 rjd 2 r± in the 
Hamiltonian expansion. 



4 Interaction of the vortex ring with its image 

Now let us for simplicity consider the case of a single ring. It is shown in the 
next section, that for a sufficiently deep ring the interaction with its image is much 
stronger than the interaction with the surface waves. So it is interesting to examine 
the motion of the ring neglecting the surface deviation. In this case we have the 
integrable Hamiltonian for the system with two degrees of freedom 

H = ^{ a ^ + p ^-^zT 1 )^ z<0 (29) 

where a ~ const. The system has integrals of motion 

P x = p = const, H = E = const 

so it is useful to consider the level lines of the energy function in the left (Z, P z )- 
half-plane taking P x as the parameter (see the Figure). 

One can distinguish three regions of qualitatively different behavior of the ring 
in that part of this half-plane where our approximation is valid (see Eq . (|20|) ) . In 
the upper region the phase trajectories come from infinitely large negative Z where 
they have a finite positive value of P z . In the process of motion P z increases. This 
behavior corresponds to the case when the ring approaches the surface. Due to 
the symmetry of the Hamiltonian (|29|) there is a symmetric lower region, where the 
vortex ring moves away from the surface. And there is the middle region, where P z 
changes the sign from negative to positive at a finite value of Z. This is the region 
of the finite motion. 

In all three cases the track of the vortex ring bends toward the surface, i.e. the 
ring is "attracted" by the surface. 
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Figure 1: The sketch of level lines of the function H(Z, P z ), Eq.(p9|). 

5 Cherenkov interaction of a vortex ring with sur- 
face waves 

When the ring is not very far from the surface and not very slow, the interaction 
with the surface waves becomes significant. Let us consider the effect of Cherenkov 
radiation of surface waves by a vortex ring which moves from the infinity to the 
surface. This case is the most definite from the viewpoint of initial conditions choice. 
We suppose that the deviation of the free surface from the horizontal plane z = is 
zero at t — > — oo, and we are interested in the asymptotic behavior of fields r\ and 
ip at large negative t. In this situation we can neglect the interaction of the ring 
with its image in comparison with the self-energy and concentrate our attention on 
interaction with surface waves only. 

The ring moves in the (x, z)-plane with an almost constant velocity. In the main 
approximation the position R of the vortex ring is given by the relations 

BE P 
R^Ct, C = C(P) = — = (C x ,0,C z )~— (30) 

C x > 0, C z > 0, t< 0. 

The equations of motion for the Fourier-components of 77 and ip follow from the 
Lagrangian (ffijp 

% = fc^k, ipk + gvk = -$k (31) 

Eliminating 77k we obtain an equation for ip^ 

4 + gktpk = -$ k (32) 
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where $k is the Fourier-transform of the function $(rj_). Simple calculations give 

$ k = / ; z xX e-« k * x+k ^dxdy = 

2ir J U x 2 + y 2 + Z 2y 



-ik x X 



2tt V z v 1 17 \Z\dk 

where 



P z D(k\Z\) + i^4-D(k\Z\)) (33) 



• V + /3 2 + l) 3 



Finally, we have for $k 

$ k = (^^^ - P 2 ^) e -*l*l-*«* = - P 2 ^) e *(*cW^) (35) 

Due to the exponential time behavior of $k(£) it is easy to obtain the expressions 
for tpk(t) and ^(t). Introducing the definition 

A k = kC z - ik x C x (36) 

we can represent the answer in the following form 

^'■ '^ ■(s)^ (37) 

= (9 5*T4 eAk ' (38) 

The radiated surface waves influence the motion of the vortex ring. The terms 
produced by the field rjk(t) in the equations of motion for the ring come from the 
part / •f]^d 2 r± in the Lagrangian (28). Using Eq.(|35|) for the Fourier-transform of $ 



we can represent these terms as follows 



SZ = f ^Li h e kz +«" x (40) 
SP. = -J ' (»*) (p. + ^) <=* W (41) 

sp ^-iw?^ k ( p ' +i ^y w <42) 



We can use Eq.(|38|) to obtain the nonconservative corrections for time derivatives 
of the ring parameters from these expressions. It is convenient to write down these 
corrections in the autonomic form 



^_( p \ f d2k ( tk A (kC z -ik x C x ) 3 2klzl 
' " 1 CJ J (2tt) 2 I k gk+ (kC z - ik x C x f 
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)• 








p 


-( 


c 



,)7 I — ) [ ^ . ( kCz lkxC *f e ~2k\Z\ (u] 

" \C)J (2tt) 2 gk+(kC z -ik x C x f V ' 

: P i - \ i ^ 2 k ( | {kC z — ik x C x f{C 2 k 2 + C%kl) 2k\z\ / 4 r\ 
■ r U * J ' (2vr) 2 I k ) gk + (kC z - ik x C x f ' ' 

(kC z — ik x C x f{C 2 k 2 + C^kp 2k\z\ /ao\ 

(2tt) 2 " gk+ (kC z -lk x C x ) 2 ' 1 ) 

where C x and C 2 can be understood as explicit functions of P defined by the de- 
pendence C(P) = dS/dP. More exact definition of C x and C z as X and Z is not 
necessary. 

To analyze the above integrals let us first perform there the integration over the 
angle ip in k-space. It is convenient to use the theory of contour integrals in the 
complex plane of variable w = cosy?. The contour 7 of integration in our case goes 
clockwise just around the cut which is from —1 to +1. We define the sign of the 
square root R{w) = y/l — w 2 so that its values are positive on the top side of the 
cut and negative on the bottom side. After introducing the quantities 

c * 2 1 . 1 [g ,._v 



we have to use the following relations 

dw w(w + ia) 3 



h(a,b) 
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\/l — w 2 b 2 — (w + ia) 2 



. (b — ia)b 2 (b + ia)b 2 , 

?r(l + 2b 2 ) + in f ' = v ; = (48) 



I 2 (a,b) = ij> —p= 



y/l - {b-iaf y/l - (-b-iaf 
dw (w + iaf 



a/1 - w 2 b 2 - (w + ia] 



2 



Ji(a, b) = i <j> 



vr 2a+ ° + " (49) 

j I -{b-iaf y/l- (-b-iaf 



dw w(w + iaf(w 2 + a 2 ) 



\/l — w 2 b 2 — (w + iaf 



-Anab 2 + 7T I b (^-^)(Q 2 + (^-^) 2 ) + b(b + ia)(a 2 + (b + iaf) . 

y/l-(b- iaf yl - {-b-iaf 



J 2 (a, 6) = y 
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(u> -Ma) (u> + a ) 
\/l — w 2 b 2 — {w + ia) 2 
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-2vr(a 2 + b 2 + 1/2) - ni 



b(a 2 + (b - ia) 2 ) b(a 2 + (b + ia) 2 ) 



l-(b-ia) 2 Jl-(-b-ta) 2 



(51) 



where the sign of the complex square root should be taken in accordance with the 
previous choice. It can easily be seen that the integrals I 2 and J\ have resonance 
structure at a < 1 and |6| < 1. This is the Cherenkov effect itself. Now the 
expressions ( |43] - f46|) take the form 



5X 



+00 



h(a,bk) k 2 e~ 2m dk 



(2vr) 2 \Cl 



FAa, 



^=(^2 / h(a,b k )k 2 e- 2k ^dk 



6P X 



p2 



(2vr) 2 \Cl 



4=\ F 2 a, 



(27T) 



Pi 



SP Z 



(2vr) 2 \Cl 



^7 Gi a 



Go [a 



Ji (a, &k) k 3 e~ 2k ^dk 

77^ / J 2 (a,6 k )^ 3 e- 2fc l z lrffc 
^ ' 

Here the functions F 1 (a, Q)..G 2 (a, Q) are defined by the integrals 

+00 . , 
F!(a,Q)= y J i^"y|j exp(-Q0^ 2 ^ 

+00 . x 
F 2 (a,Q)= | J 2 ^a,- 7 =Jexp(-gOe 2 ^ 

Gi(a,Q) 



2gg| ' 

eg , 
2^1 

2g\Z\ 
C 2 X 



(52) 

(53) 
(54) 
(55) 





+00 





+00 



J i(«,-y|] exp(-QOe 3 ^ 



G 2 (a,Q) 



Ji a, 



exp(-QO£ 3 d£ 



(56) 

(57) 
(58) 
(59) 



and Q = 2g\Z\/C 2 is a dimensionless quantity Q. The Cherenkov effect is most clear 
when the motion of the ring is almost horizontal. In this case a — > +0, and it is 
convenient to rewrite these integrals without use of complex functions 

+00 1 
fi(+0, Q) = vr I (^ 2 + 2£) exp (-Q0 d£ - 2vr j -^=L= exp (-QO (60) 



3 If we consider a fluid with surface tension a, then two parameters appear: Q and T = ga/C^- 
In that case one should substitute — > + T£ as the second argument of the functions 

Ii,h,Ji, J2 in the integrals (|56- 59 ) 
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F 2 (+0, Q) = G^+O, Q) = -27T j exp (-Qf) (61) 



G 2 (+0, Q) = -7T / (e 3 + 2£ 2 ) exp (-gO ^ + 2tt y -^=J= exp (-QO (62) 



Here the square root is the usual positive defined real function. We see that only 
resonant wave- numbers contribute to the functions F 2 and G±, while F 1 and G 2 are 
determined also by small values of £ which correspond to the large scale surface 
deviation co-moving with the ring. So the effect of the Cherenkov radiation on the 
vortex ring motion is the most distinct in the equations for Z and P x . Especially 
it is important for P x because the radiation of surface waves is the only reason for 
change of this quantity in the frame of our approximation. 

The typical values of Q are large in practical situations. In this limit asymptotic 
values of the integrals above are 

Fl(+ °' g) ^"2^' <? 2 (+0,Q)«^ 

exp(-Q) 



F 2 (+0,Q) = Gi(+0,Q) « -2%^ 



VQ 



and 



1 p2 q p2 I 

SP X « =— ■ Q 3+1/2 exp(-Q), SP Z w + — ^r • -. 

32v^F|Z| 4 v PV 32tt|Z| 4 Q 

It follows from these expressions that the interaction with the surface waves is small 
in comparison with the interaction between ring and its image, if Q ^> 1. The 
corresponding small factors are 1/Q for X and P z , and Q 2+1 ^ 2 exp(— Q) for Z. As 
against the flat boundary, now P x is not conserved. It decreases exponentially slowly 
and this is the main effect of Cherenkov radiation. 

We see also that the interaction with waves turns the vector P towards the surface 
which results in a more fast boundary approach by the ring track. 



6 Conclusions and acknowledgments 

In this paper we have derived the simplified Lagrangian for the description of the 
motion of deep vortex rings under free surface of perfect fluid. We have analyzed the 
integrable dynamics corresponding to the pure interaction of the single point vortex 
ring with its image. It was found that there are three types of qualitatively different 
behaviour of the ring. The interaction of the ring with the surface has an attractive 
character in all three regimes. The Fourier-components of radiated Cherenkov waves 
were calculated for the case when the vortex ring comes from infinity and has both 
horizontal and vertical components of the velocity. The non-conservative corrections 
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to the equations of motion of the ring, due to Cherenkov radiation, were derived. 
Due to these corrections the track of the ring bends towards the surface faster then 
in the case of flat surface. For simplicity, all calculations in Sec. 5 were performed for 
a single ring. The generalization for the case of many rings is straightforward. 

The author thanks professor J.J. Rasmussen for his attention to this work and for 
helpful suggestions. This work was supported by the INTAS (grant No. 96-0413), 
the Russian Foundation for Basic Research (grant No. 97-01-00093), and the Landau 
Postdoc Scholarship (KFA, Forschungszentrum, Juelich, Germany). 
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